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Abstract. On the basis of the Kadanoff-Baym (KB) version of the time-dependent Green's 
function method, a new Ansatz for the approximation of a spectral function is offered. The 
Ansatz possesses all the advantages of quasiparticle and extended quasiparticle approxima- 
tions and satisfies the KB equation for a spectral function in the case of slightly nonequi- 
librium system when disturbances in space and time are taken into consideration in the 
gradient approximation. This feature opens opportunities for the microscopic derivation of 
the Landau kinetic equation for the quasiparticle distribution function of the normal Fermi 
liquid and provides the widening of these equation's temperature range of validity. 
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1. Introduction 

The only microscopic theory which is capable of describing both the statistics and dynam- 
ics in a comprehensive way is the Green's function approach developed in different varieties 
for equilibrium and nonequilibrium problems and for zero and finite temperatures. The ini- 
tial emphasis in application of Green's function techniques was to understand the properties 
of normal condensed matter systems, superconductors, and superfluids. 

Microscopic models used for a description of quantum interacting many-body systems in- 
volve spectral functions which play a central role in the formalism. Spectral functions are 
fundamental in describing the nuclear correlations, electrons correlations in metals and semi- 
conductors, and in many other equilibrium and nonequilibrium properties. The most general 
approach to the problem is based on the real-time Green's function formalism of Martin and 
Schwinger, further developed by Kadanoff and Baym [lj. Keldysh j2] developed a diagram 
version of the theory for nonequilibrium systems equivalent to the Kadanoff-Baym (KB) 
approach. The KB equations were used in several contributions to this evolving field with 
applications to nuclear matter [3J, to one- and two-band semiconductors [H E], to electron 
plasmas [6] [E], etc. The various approximations of the KB equations differ essentially by 
the reduction schemes of the two-time Green's functions to the reduced equilibrium density 
matrix and to the quasiparticle distributions. There are several papers on spectral functions 
published during the past years using different methods and approximations [9] [13] . The 
nonequilibrium extension of the KB formalism has been recovered within the quasiparticle 
approach to kinetic equation for weakly interacting particles and referred to as a modified 
KB Ansatz [22] [2H - 

In a variety of works beginning with pQ, the spectral function was approximated by a 
delta function of energy peaked at the quasiparticle energy. Such a simple approximation 
falls short of describing many important features of the systems under consideration. Im- 
proved spectral functions were necessary and some variants were offered in [TT1 [121 [23]. The 
extended quasiparticle (EQP) approximation was introduced and it was shown to remedy 
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some of the faults of the simpler quasiparticle (QP) approximation . Thus the EQP approx- 
imation was a considerable improvement of the formalism and it was used to compare the 
zero-temperature case of the KB approach with Brueckner theory of nuclear matter (TTJ [12] . 
The EQP approximation can be modified slightly while maintaining its simplicity and renor- 
malization property. One possible modification of EQP was offered in |12j . but it was not 
even mentioned among possible forms of spectral functions in [13j . 

Further development of the KB theory was performed in several different directions. Nu- 
merical solutions of the KB equations were presented in [H]. The extended quasiparticle 
picture which was first offered in |15j . was developed for small scattering rates in [T6]~|18]. 
The separation of pole and off-pole parts of a spectral function was discussed in detail in [TF] . 
The nonlocal quasiparticle kinetic equation for the momentum-, space-, and time-dependent 
distribution function, which has the form of a Boltzmann equation with the quasiparticle 
energy, was derived in [18]. The generalization of the KB equations to the case of arbitrary 
initial correlation was performed in [19] . Correlation effects related to the collision integral 
were shown to cause a damping of the spectral function [2U]. However, for the present analy- 
sis (see below), like in some other cases [2T], a detailed self-consistent microscopic treatment 
of the correlation effects on the spectral function can be avoided. 

Besides other applications, the spectral function was used for the derivation of the kinetic 
equation of the phenomenological Landau normal Fermi liquid theory. The initial derivation 
of this equation was produced by Kadanoff and Baym on the basis of the quasiparticle 
approximation for the spectral function and was continued by some followers who used the 
extended quasiparticle approximation. In both cases, the second Poisson bracket in the right 
side of the KB generalized kinetic equation could not be eliminated in a lawful mathematical 
way. This was the reason of the narrowing of the temperature range of validity of Landau's 
equation. 

This paper is devoted to the explanation of the successfulness of the phenomenological 
Landau theory of normal Fermi liquid and its applicability far beyond its temperature range 
of validity established on the basis of quasiparticle and extended quasiparticle Ansatze for 
the spectral function. We offer another approximation for the spectral function which pos- 
sesses the advantage typical for QP and EQP approximations but has an additional property 
to satisfy exactly the KB equation for the spectral function in the case of slightly nonequi- 
librium systems. It makes this form of the spectral function preferable in the nonequilibrium 
systems and opens certain opportunities for the widening of the range of validity of Landau 
Fermi-liquid kinetic equation. 

The following section contains a concise presentation of the necessary formulas of the KB 
theory for both equilibrium and slightly nonequilibrium systems which will be used in Sec. 
3 for a comparison of three different approximations for the spectral function. Section 4 
contains the analysis of the question that what approximation for a spectral functions suits 
better to the strict results of Kadanoff and Baym for the nonequilibrium case and deals with 
the Landau kinetic equation for the quasiparticles in the normal Fermi liquid which is proved 
to be valid in a wider temperature region than it was considered in the initial derivation on 
the basis of the KB formalism. In Sec. 5 a discussion and summary are presented. 



2. Main formulas of the KB formalism 

We will introduce the main results of the KB theory for fermion systems [1] keeping in 
mind the questions which will be discussed below. The KB formalism leads to the following 
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general expression for the one-particle spectral function a(pu) of a system in equilibrium: 

T(pu) 



(1) a(pu) 



[u - E HF {p) - Re a c (pu)} 2 + 



T 2 {pu}) ' 



where E HF (p) is a one-particle energy in the Hartree-Fock approximation. Real and imag- 
inary (r) parts of the correlation self-energy function a c are related to each other through 
the Hilbert transform, 



(2) Re a c (poo) = P J 



duJ T(pu) 



oo 



2lT U — u' 



Here, P refers to a principal value integration. 
The spectral function satisfies the exact sum rule: 

f°° du N 

(3) y — o(^) = i, 

for all the values of p. This result follows directly from the commutator relations for field 
operators and can serve as a keystone for the checking of the validity of all approximations 
for the spectral function (1). 

In the Hartree-Fock approximation, when the correlation self-energy function Re a c (pu) = 
0, the spectral function turns to be a delta function of the Hartree-Fock energy: 

(4) a HF (pu) = 2tt8[u - E HF {p)l 

and the sum rule (3) is trivially satisfied identically for any value of a HF (p). The simplest 
approximation in the case when Re a c (puj) ^ is the so called (QP) approximation: 

(5) a QP = 2TrZ(p)5(u-E(p)), 
where E(p) is the solution of the equation: 

(6) E(p) = E HF (p) + Re a c [p, E(p)}, 
and the renormalizing factor Z(p) is given by the expression 



(7) z-\p) = i - 

In the QP approximation, the sum rule (3) reads 



u=E(p) 



(8) J ^-a QP (puj) = Z(p). 

A severe drawback with oqp is that it normalizes to Z(p) rather than to 1 as in Eq. (3). We 
will see below how this deficiency is removed by the extended quasiparticle approximation 
[TTJ [13] and by other more advanced approximations. 

In the case of slowly varying in space and time disturbances, after the transition to Wigner 
coordinates, 



(9) R= -(rl + ry), r = rl-ry , T=~{t 1 + t 1 >), t = t 1 -t 1 >, 



-(ri + Tv), r = n-r 1 / , T= -( 

and the performance of the Fourier transform with respect to r, t, all the quantities entering 
the theory are considered to be the functions of p, to, R, T, for example, 



a = a{pu>\ RT). 
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If we take into account only the first derivatives with respect to slowly varying quantities R 
and T in the KB equations for the correlation functions, we come to the following equation 
for the spectral function a(pu; RT) [Tj: 

[u - E HF {p; RT) - Re a c {pu; RT),a{puj; RT)} 

(10) + [Re g(pco; RT), T(pu; RT)} = 0, 

and to the generalized KB kinetic equation for the correlation function g < (pu; RT): 

[u - E HF (p; RT) - Re a c (pu; RT), g < (pu; RT)} 

(11) + [Re g(poo;RT),a < (pu;RT)} = (a < g > — a > g < )(puo; RT). 

Here, [A, B] is the so called generalized Poisson bracket, defined by the expression: 

< 12 > l^'^f-if -v^.v, B + v^ .V ?B , 

E HF and Re a c include the interaction with the external field U(RT). The exact solution of 
Eq. (10) is given by the expression 

(13) g(pz; RT) = [z - E HF (p; RT) - Re a c (pz; RT)}' 1 , 

where all the functions, such as a(pu; RT), Re g{pu; RT), etc., are determined by the same 
formulas as in the equilibrium case with all entering quantities depending on p, u, R, T. For 
example, 

J — C 



-oo 



2-7T UJ — UJ 



In fact, the solution (13) gives almost the same evaluation of the spectral function a as in 
the equilibrium case: 

a(pu;RT)- 



[u - E HF (p; RT) - Re a c (pu; RT)} 2 + E^E^Il 

Since the solution (13) is of exactly the same form as the equilibrium solution, it must 
reduce to the equilibrium solution as T — > — oo. Thus, it satisfies the initial condition on 
the equation of motion. This result means that the approximation for the nonequilibrium 
spectral function can be written in the same form as in the equilibrium case. 

Equation (11) provides an exact description of the response to slowly varying disturbance. 
All the quantities appearing in this equation may be expressed in terms of the correlation 
functions g < and g > . In particular, a < and a > are defined by a Green's function approxi- 
mation that gives the self-energy in terms of g K and g > . 

3. Approximations for the spectral function 

We will consider three different approximations for the spectral function which we will 
call Ansatz 1, Ansatz 2, and Ansatz 3 correspondingly. Thus, we have 

(14) Ansatz 1 : a EQP = 2ixZ{p)b{uj - E) + ; 

(15) Ansatz 2 : a iQP = 2nZ{p)5{uj - E) + Z{p)P^^; 

(u-E) 2 ■ 



(16) Ansatz 3 : a { = 2nZ{p)5{uj - E) + Z 2 {p)P 1 ^- 



Here, E is determined by Eq. (6). 
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The Ansatz 1 corresponds to the extended quasiparticle approximation which was in- 
troduced, discussed, and used for numerical calculation in pj] p3]. The first form of this 
Ansatz was offered in [25J. The Ansatz 2 as it was already mentioned above, offered in 
[12] . but not discussed and was never used for the numerical calculation and was not even 
mentioned in [T3]. A little bit different in technical details but the same in principle way 
of introducing the Asatz 1 was offered in [23J. We would like to emphasize that the forms 

(14) and (15) of the spectral function were constructed in [H] [H], [23J on the basis of the 
general expression (1) and were obtained by means of the Taylor expansion in powers of T 
in the frame of different approximations. In reality, there does not exist a mathematically 
strict correct form for the expansion of (1) in power series of r(o>) which starts with the 
delta function when T(u) = 0. Thus, expression (14) and (15) should be considered as some 
true-like approximate forms, as it is recognized in [23J , which have better qualities compared 
with the quasiparticle Ansatz (5). In particular, Eqs. (14) and (15) obey the sum rule (3) 
(see below). Neither can be obtained in a strict way the improved Ansatz 3 dj offered by 
us. The origin of the possibility of such Ansatz on equal terms with expressions (14) and 

(15) can be shown on the basis of the following consideration. 

We start with a well known relation of the Fourier transform in the case of a constant 
value of c: 

(17) / e~^ c e ltx dt = - , c> 0. 



C_ I ry*£. 
r J/ 

Expanding the first exponent in the left side of Eq. (17) in Taylor series, we get 

2r f°° / r 2 t 2 \ 



Now we use the formulas equivalent to those represented in 

/■OO 

(19) / t {2n) e ttx dt = 2n(-i) 2n 5 (2n) (x), n = 0,1,2, 



(20) 



/OO , . -| 

\t\ {2n+1) e ltx dt = -2sin((2n + l)-)(2n + l)! — 
-oo v 2/ M 



2n+2 ' 



n = 0,1,2,... . 



If the quantity T(u) in Eq. (1) is constant, then the expression (18) would lead to a strict 
correct expansion of spectral function a(u) in terms of the power series of T. However, T(u) 
cannot be constant due to the dispersion relation (2). In the case of T(u) depending on u, 
one can rely only on the first two terms of the expansion: the delta function independent of 
T(u) and the term proportional to T(u). Then one should notice that the quantity x in Eq. 
(17) or (18) corresponds tow - E HF (p) — Re a c (puj) in the formula (1). Taking into account 
the relation (6), it becomes clear with the precision to the second derivative 9 ^ 2 ° c , the 
term of the expansion (18) proportional to T(u) involves Z 2 (p) in the numerator, and we 
come to the formula (16). Indeed, in the mentioned approximation we have 

u - E HF (p) - Re a c (puj) = Z~ x {$)[u - E(p)]. 

Now we should compare the expressions (14)-(16) on the basis of strict results obtained in 
the KB theory. We will check the correspondence of these expressions to the sum rule (3), 
the second sum rule on the basis of energy considerations and to Eq. (10) which determines 
the spectral function in the case of slowly varying in space and time disturbances. We start 
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with the sum rule (3). Substitute Ansatz 1 given by Eq. (14) to Eq. (3) and we get 

duo T(u) 



2tt (to - Ef 



(21) jH ^a EQP = Z{p) + Pj 

Due to the dispersion relation (2) we have 
(00 , d Re a c (u) , r <U< L-(-' 



du '* 7„ 0O 2vr(a;-i?)2' 

and Eq. (21) becomes 

y°°da; 1 9Rea c H 

(23) / —a E Qp 



DO 



2vr ^ i- ^^M | F 5w 



with the precision to v duj 
The Ansatz 2 gives 



9Reac\2 



—a lQP = Z(p)(l + P J 



du „ r 00 rfa; r(u) 

27r (a; - £^) ; 

Taking into account the relation (22), we see that in the case under consideration the integral 

(24) is equal to 1. 

Substituting Eq. (16) into Eq. (3) and taking into account (22), one gets 

(25) y_ 0i = Z( g)(i_ Z _ii !,)„!, 

with the precision to ( 9 Rc 9 ° c ^ ) 2 - 

Thus, formally the Ansatz 2 is the best in the sense of sum rule (3), but errors brought 
by Ansatz 1 and Ansatz 3 can be neglected in the approximation under consideration. 

We can check the expressions (14)-(16) also with the help of the second sum rule 

(26) I" pua(pu) = E HF {p). 



oo 



2tt 



The first terms in the right side of Eqs. (14)-(16) give the same result 

(27) r —u2nZ(p)5(u - E) = Z[E HF + Re a c (E)]. 

J-oo 2vr 

The contributions of the second terms in right sides of Eqs. (14)-(16) differ only in factors 
depending on the power of the renormalizing factor Z. The corresponding integral in each 
case is: 

du; uTjuj) /-°° du I» Ep r dco T(u) 



2tt (co - Ef 2n (oj - E) J _ x 2tt (u - Ef 

(28) = Z-\\ - Z)E HF -{2Z -l)Re a c {E)]. 

The last equality in Eq. (28) follows with the help of the expressions (2), (7), and (22). 
Now taking into account the renormalizing factors standing with the second terms in Eqs. 
(14)-(16), we get with the help of Eqs. (27) and (28) that the second sum rule (26) is 
valid for all the ansdtze (14)-(16) with the same precision up to the terms of the order 
( - R g^ c ^ ) 2 and in this sense they are equivalent in principle. However, they turn to be not 
equivalent in the sense of satisfying the Eqs. (10) and (11), although it is mentioned in [23] 
that the corresponding off-pole part of g K for the Ansatz 1 (14) compensates a "dominant 
part of the puzzling term [Re 5f,<r < ]" in the generalized kinetic equation (11). However, 
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represented in [23] analysis did not lead to mathematically lawful complete elimination of 
this puzzling term. In the next section, we show that the Ansatz 3 (16) solves the problem 
and consequently is the best in the nonequilibrium case. 

4. Kinetic equation of the Landau theory of the Normal Fermi Liquid 

According to the strict result (13), every approximate Ansatz for the spectral function 
should keep valid in the case of slowly varying in space and time disturbances. The only 
difference occurs due to the dependence of all quantities on R and T. When we substitute 
Eqs. (14)-(16) to Eq. (10), the first terms in the right side of these expressions give 

[u - E HF (p) - Re a c (pu; RT),2irZ(p)5(uj - E)\ 



[oo 



E) 1- 



d Re cr r 



du 



5{uo - E) 



(29) = 2nZ(p) 

due to the property of the generalized Poisson bracket 

(30) [A, f(A)} = 0. 

The second term in the right side of Eq. (14) leads to the expression 

r(jxo;RT)- 



0. 



w - E Ht (p) - Re a c (pco; RT), 



[CO 



Ef 



(31) 



z-\p)- 



CO 



E,T], 



(lu-E) 21 

with the precision to - K qJi c ^ ■ The second Poisson bracket in Eq. (10) due to Eq. (13) 
gives with the same precision the expression 

[Re g{pu- RT) , T(pu; RT)} = Z(p) f ' 



[CO 



(32) 



-m- 



UJ 



-e,t\. 



[to - E) 2 

Finally, collecting all the terms, we get that the left side of Eq. (10) in the case under 
consideration is different for Eqs. (14)-(16) and equals to 



(33) 
(34) 



Ansatz 1 



Ansatz 2 



Ansatz 3 



E) 



\-\LO 



E) 



r 



[Z(p) - Z-\p)); 

[zip) - 1]; 



[zip) - Z{p)\ = o. 



.ico-EY 

The Ansatz 3 satisfies Eq. (10) exactly. It means that this Ansatz is preferable for con- 
sidering slightly nonequilibrium systems when disturbances slowly vary in space and time. 
In particular, it turns out that the Ansatz 3 (16) opens opportunities for the derivation of 
the Landau Fermi liquid's theory kinetic equation for the quasiparticle distribution function. 
We will see that the usual considerations about the smallness of the quantity T(o;) are not 
necessary anymore and the puzzling term [Re 5',cr < ] in Eq. (11) is canceled completely by 
the off-pole part of the approximation (16). 

From the pioneer work of Kadanoff and Baym, the derivation of the kinetic equation of 
normal Fermi liquid theory was based on the assumption of the smallness of the functions 
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a < and a > near the Fermi level \x of the system [U [29] . Indeed, defining a local occupation 
number f(puj; RT) by writing 

(35) g<{pu- RT) = a(pu; RT)f(puj; RT), 
where in equilibrium, at zero temperature 

(36) /(^;*TWM={; % 
one comes to the relations 

cr > (pLj) =0 for uj < ji , 

(37) a<(pcj)=0 for w > // . 

The proof of these relations depends only the fact that / = 1 for uj < p, and / = for 
uj > \i. Since we are assuming that / has a similar behavior in the slightly nonequilibrium 
case, it follows that 

a > (puj; RT) = for uj < li(RT) , 

(38) a < {puj; RT) = for uj > fJ,(RT) , 

where li(RT) is a local chemical potential of the system. For normal Fermi systems, a > and 
cr < are continuous functions at oj = [i. It means that r = cr > + a < is small near il. As it was 
shown in [1], these assumptions about / for uj appreciably greater or less than n(RT) lead 
to a consistent solution of the generalized kinetic equation (11). In such a situation when 
a > and cr < are both negligible for u near fx, Eq. (11) in collisionless approximation can be 
written in the form 

(39) [uj - E HF (p; RT) - Re a c (p; RT),a(puj; RT)f(pu; RT)] = 0. 

The second Poisson bracket in the right side of Eq. (11) is dropped out due to the discussed 
considerations about the smallness of F near /j,. 

Using the quasiparticle Ansatz (5) for the spectral function a, one comes to the kinetic 
equation for the quasiparticle distribution function n(p, RT) of the normal Fermi liquid 
theory pQ: 

(40) — + V P E. V n n-V M E. = 0, 

n(p; RT) = f(puj; RT) | u=BtfljAr) • 

The smallness of the function T is a keystone in the above considerations and such a smallness 
looks like a necessary condition of the validity of Eq. (40), which is supposed to be valid only 
in the region uj w fi(RT). However, it can be shown that neglecting collisions, Eq. (40) stays 
valid with the precision to T 2 if one uses Ansatz 3 (16) for the spectral function a instead 
of the quasiparticle Ansatz (5). In this case, the collisionless equation (11) is written in the 
form 

(41) [uj - E HF - Re a c , af] + [Re g, a<] = 0. 

The first term in all the Ansatze (14)-(16) which coincides with the quasiparticle Ansatz 
(5) leads directly to Eq. (40) as it is shown in pp. The second Poisson bracket in Eq. (41) 
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in the linear approximation in T gives the result 

1 



[Re g, o < 



- E HF - Re a r 



,17 



z 



uj-E 



,17 



(42) 



Z[Tf,u-E]- 



u - Ef 



It is easy to see that this term will be compensated by the second term in the Ansatz 3. 
Indeed, substituting this second term to the first Poisson bracket in Eq. (41), one gets 



(43) 



uj-E 



HF 



Re a r 



17 



{u - E)< 



Z[u-E,Tf] 



(uj-E) 2 ' 



and this term compensates Eq. (42) in Eq. (41). Observe that neither Ansatz 1 nor Ansatz 
2 will lead to this result. 

We would like to stress that all considerations based on the relations (36)-(38) were 
necessary only for the sake of the elimination of the second Poisson bracket in Eq. (41). The 
usage of the Ansatz 3(16) makes the smallness of the function T not necessary. The quantity 
T(uj) can be finite. For qualitative estimation of the precision, we can use the third term in 
the expansion (18) which corresponds to T 2 (uj). Substituting the term proportional to T 2 (uj) 
into Eqs. (10) and (11), it is not difficult to show that in the gradient approximation these 
equations are valid up to the terms of the order r 2 (u;). 

This result means that the collective excitation spectrum of the system which is determined 
by the kinetic equation (40) is not changed when we take into considerations the term in 
the spectral function proportional to T. At the same time, physical quantities which are 
determined by the correlation function g K will obtain additional terms due to the second 
term in the expression for the spectral function. In particular, the reduced density matrix, 

(44) p(p; RT) = J°° |V &>\ RT) , 

will possess power tails while the quasiparticle distribution n decrease exponentially with 
energy. 



5. Discussion and Summary 

The differences in the form of the expressions (14)-(16) come through the ways they were 
obtained. When the Ansatz for the spectral function is obtained by a Taylor expansion of 
Eq. (1) around the quasiparticle peak, it involves the necessity of the conditions T << Re a c 
and 9 g^ CTc << 1 [H]-[l3]. The approach (18) based on the Fourier transform demands 

another condition (u — E) a << 1, which means that a relation between T and Re a c 

can be arbitrary, in the frame of the Hilbert transform (2). 

The possibility for V to be not small plays an important role for the increasing of the range 
of validity of the kinetic Fermi liquid equation. The comparison of the results obtained on 
the basis of this equation with experimental data confirmed the validity of the equation far 
beyond the limits established by the conditions (36)-(38) and, consequently, for a essentially 
larger temperature interval. 

The Ansatz 3 (16) leads to the elimination of the second Poisson bracket in Eq. (11) in 
a very natural way. As it turned out, this term in the generalized KB kinetic equation (11) 
is the main obstacle on the way of the extension of Fermi liquid equation (40) to a larger 
temperature interval [23]. We should mention that the validity of the kinetic equation (40) 
in the case of finite values of V was shown by different mathematical method in [30J but the 
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result looked contradicting to the commonly accepted conditions (36)-(38). Only the direct 
"lawful" elimination of the term [Re g, <r < ] in Eq. (11) makes the situation clear. 

Despite the absence of the rigorous foundation of the validity of the Fermi-liquid equations 
in the case of large T, these equations were successfully used for the variety of systems of 
strongly interacting particles. Thus, such equations were used for theoretical description of 
the magneto-ordered state in 3d metals in the framework of the Anderson periodic model 
|31j . The value of T was considered to be constant and of the same order as other energy 
parameters of the system: the width of the d level was considered to be equal to a relatively 
large jump parameter v : v ~ V ~ 1 eV. The obtained results turned to be in a good 
agreement with experimented data for the photoeffect on polarized electrons [32, 33J. In 
conclusion, we will mention that the Ansatz 3 (16) can be preferable also for calculations of 
the equilibrium properties such as produced for a nuclear matter in |12j . As it is stated in 
this work, using the extended quasiparticle approximation (14), the total binding energy is 
obtained to be 

E EQP = 17 A MeV/nucleon. 
In contrast, the quasiparticle approximation (5) leads to the result 

E QP = 15.9 MeV/nucleon. 
Bruckner theory meanwhile gives 

Eg = 16.7 MeV/nucleon. 

Due to the extra factor Z 2 (p)[Z(p) < 1] in the second term of the Ansatz 3 compared with 
the Ansatz 1, it is clear that the result of this Ansatz must be between the values of Eqp 
and Eeqp- It will be closer to the value of Eb than Eeqp and closer to the value of Ep than 
Eqp- 
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